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Abstract 

The  GPS  system  is  undergoing  a  modernization  effort  which  will  add  several 
new  signals  to  be  placed  at  the  LI,  L2,  and  L5  frequency  bands.  One  of  the  signals 
to  be  placed  on  L2  is  a  new  military  code  (M-code)  which  may  be  transmitted  at  a 
higher  power  level  than  current  GPS  signals.  Other  users  of  the  L2  frequency  band 
are  concerned  with  the  potential  interference  that  may  be  caused  by  the  increase  in 
power  of  the  GPS  signal.  One  particular  use  of  the  1215-1400  MHz  frequency  band 
is  Air  Traffic  Control  (ATC)  radar  applications. 

This  thesis  models  an  ARSR-4  radar  system  as  a  single  pulse  radar  and  simu¬ 
lates  the  effects  that  the  M-code  signal  has  on  a  radar  system’s  detection  capability. 
Looking  at  a  worst  case  radar  detection  scenario,  where  the  M-code  signal  power 
incident  on  the  radar  would  be  at  a  maximum,  the  results  indicate  that  for  the  min¬ 
imum  specified  received  M-code  signal  power,  the  effects  are  minimal.  However,  for 
the  maximum  specified  M-code  signal  power,  the  effects  are  quite  noticeable  to  the 
point  where  the  M-code  may  prevent  detection  of  a  single  pulse. 

The  results  only  apply  to  a  worst-case  scenario  for  a  single  pulse  radar.  Further 
modeling  of  the  ARSR-4  system,  to  include  pulse  integration,  would  be  necessary 
before  a  definitive  conclusion  can  be  drawn  that  the  M-code  will  significantly  affect 
radar  system  performance. 
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THE  EFFECTS  OF  GPS  M-CODE  ON  RADAR  DETECTION 


I.  Introduction 

In  a  letter  dated  10  August  2001,  the  International  Civil  Aviation  Organization 
(ICAO)  raised  the  following  issue  to  be  addressed  at  the  International  Telecommu¬ 
nication  Union  (ITU)  World  Radiocommunication  Conference  (2003)  (WRC-2003): 

“...the  need  to  secure  the  protection  of  distance  measuring  equipment 
(DME)  and  radar  systems  from  interference  from  radio  navigation  satel¬ 
lite  service  (RNSS)  in  the  frequency  bands  960-1215  MHz  and  1215-1400 
MHz”  [16] 

The  ICAO  is  concerned  about  the  potential  interference  that  may  result  from 
the  GPS  modernization  effort.  The  current  received  GPS  signal  power  is  below  the 
noise  floor  so  that  most  systems  receiving  the  GPS  signal  as  interference  may  not  be 
able  to  distinguish  the  actual  signal  from  noise.  However,  the  modernization  plan 
calls  for  an  increase  in  the  received  GPS  power  of  some  signals  by  levels  20  dB  greater 
than  the  current  signals  [14].  The  increase  in  power  pushes  the  received  GPS  signal 
close  to  the  noise  floor,  which  may  cause  interference  to  systems  operating  near  the 
GPS  frequency.  This  concern  is  illustrated  in  Figure  1.1,  which  shows  the  proposed 
GPS  military  signal  (M-code)  at  two  received  power  levels  in  relation  to  the  noise 
floor. 

1.1  Background 

The  Global  Positioning  System  (GPS)  was  designed  in  the  early  1970s  to  pro¬ 
vide  military  and  civil  users  with  position  estimates  using  radio  navigation  (radio 
waves)  from  a  constellation  of  satellites  orbiting  20,000  km  above  the  earth.  Since 
then,  the  world  has  changed,  but  the  GPS  policies  and  regulations  adopted  by  the 
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Figure  1.1  Received  M-Code  Power  in  Relation  to  the  Noise  Floor 

U.S.  Government  for  military  and  civilian  users  have  generally  remained  constant. 
An  increasing  emphasis  on  commerce  and  the  world’s  current  political  climate  has 
given  the  civil  users  and  GPS  industry  in  the  United  States  new  influence  to  push 
for  modifications  in  GPS  design  and  policies.  This  emphasis  became  evident  with 
the  deactivation  of  Selective  Availability  (SA)  in  May  2000,  which  marked  the  first 
important  change  in  U.S.  policy  [7]. 

The  deactivation  of  SA  was  the  result  of  the  March  1996  Presidential  Decision 
Directive  (PDD)  [3]  which  announced  the  “intention  to  discontinue  the  use  of  GPS 
Selective  Availability  (SA)  within  a  decade”.  The  primary  purpose  of  SA  was  to 
intentionally  degrade  GPS  to  reduce  the  capability  of  enemy  forces  using  GPS.  To 
replace  the  function  of  SA,  the  PDD  also  tasked  the  Department  of  Defense  (DoD) 
“to  develop  measures  to  prevent  the  hostile  use  of  GPS  and  its  augmentations  to 
ensure  that  the  United  States  retains  a  military  advantage  without  unduly  disrupting 
or  degrading  civilian  use.” 
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To  comply  with  the  directive,  the  GPS  Joint  Program  Office  (JPO)  and  the 
Space  Warfare  Center  (SWC)  conducted  a  study  to  identify  the  best  GPS  signal 
redesign  to  meet  the  needs  of  both  the  military  and  civilian  GPS  users  for  the  next 
30  years  [2],  The  results  of  that  study  were  incorporated  into  what  is  now  termed 
GPS  Modernization. 

1.1.1  Current  Signal  Structure.  The  current  GPS  signals  reside  on  two 
L-band  frequencies:  LI  at  1575.42  MHz  and  L2  at  1227.6  MHz.  The  LI  frequency 
consists  of  two-spread-spectrum  BiPhase  Shift  Keyed  (BPSK)  modulated  signals  in 
phase  quadrature:  the  Precise  (P)-code,  with  a  10.23  MHz  chipping  rate,  and  the 
Coarse/ Acquisition  (C/A)-code,  with  a  1.023  MHz  chipping  rate  [2],  The  P-code  can 
be  (and  currently  is)  encrypted  for  anti-spoofing.  When  encrypted,  the  P-code  is 
denoted  P(Y)-code.  The  C/A-code  is  unencrypted  and  is  used  by  military  receivers 
for  initial  acquisition,  while  for  civilian  receivers  it  is  the  only  navigation  signal.  The 
L2  frequency  currently  consists  of  only  the  P(Y)-code. 

1.1.2  Modernized  Signal  Structure.  The  modernization  plan  calls  for 
reusing  the  existing  frequency  bands.  The  L-band  frequencies  are  especially  well- 
suited  for  radionavigation  from  space  and  additional  L-band  spectrum  for  GPS  is 
scarce.  In  addition,  the  technical  limitations  and  the  costs  involved  in  modifying 
satellites  to  transmit  at  frequencies  different  from  LI  and  L2  are  prohibitive  [2], 
Therefore,  the  modernized  signal  architecture  will  include  the  current  signal  struc¬ 
ture  with  the  addition  of  two  new  civil  signals— one  on  L5  centered  at  1176.45  MHz 
and  another  on  L2  (L2C).  In  addition,  a  new  military  signal  (M-code)  will  be  placed 
on  both  LI  and  L2.  The  resulting  modernized  GPS  signal  architecture  is  shown  in 
Figure  1.2. 
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New  Civil  Signal  C/A  Code  Signal  Y  Code  Signal 
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Figure  1.2  Modernized  GPS  Signal  Architecture  [4] 

1.2  Problem  Statement 

Since  GPS  Modernization  includes  two  new  signals  on  the  L2  band,  other 
users  of  this  frequency  spectrum  are  concerned  with  the  potential  Radio  Frequency 
Interference  (RFI)  that  may  be  caused  by  the  modernized  GPS  signals.  To  allay 
any  concerns  with  the  modernization  plan,  the  interference  issue  must  be  addressed. 
Therefore,  the  objectives  of  this  thesis  are  to  model  and  simulate  the  effects  of  the 
modernized  GPS  signals  on  radar  systems  operating  in  the  L2  band. 

1.3  Summary  of  Current  Knowledge 

Numerous  studies  have  examined  the  effects  of  RFI  on  GPS.  Most  conclude 
that  GPS  is  vulnerable  to  both  intentional  (jamming)  and  unintentional  interference 
sources.  This  vulnerability  is  a  byproduct  of  the  Direct  Sequence  Spread  Spectrum 
(DS/SS)  GPS  signal  structure  and  the  long  baseline  between  the  satellites  and  re¬ 
ceivers,  which  results  in  extremely  low  received  power  levels. 

Very  little  attention  has  been  paid  to  the  effects  of  GPS  interfering  with  other 
applications.  As  previously  stated,  the  current  GPS  signal  power  is  below  the  noise 
floor  and  most  systems  receiving  the  GPS  signal  as  an  interferer  may  not  be  able  to 
distinguish  the  actual  signal  from  noise.  Furthermore,  the  National  Telecommunica¬ 
tions  and  Information  Administration  (NTIA)  has  stated  that  at  the  current  GPS 
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signal  power  level,  which  has  not  changed  since  1978,  there  have  been  no  reports 
of  harmful  interference  to  other  users  in  the  frequency  bands  used  by  GPS  [13]. 
However,  with  the  increase  in  M-code  power,  this  situation  may  no  longer  be  the 
case. 

The  Federal  Aviation  Administration  (FAA)  recently  conducted  tests  to  de¬ 
termine  what  affect  an  increase  in  GPS  satellite  power  would  have  on  some  of  their 
radar  systems.  By  varying  the  Interference-to-Noise  (I/N)  ratios  from  approximately 
-18  dB  I/N  to  +21  dB  I/N,  the  tests  were  designed  to  determine  the  impact  on  prob¬ 
ability  of  detection,  reduction  in  range,  and  the  desensitizing  of  the  radar  [15].  As 
of  Feb  2003,  the  final  FAA  test  report  has  yet  to  be  completed. 

1-4  Scope 

As  discussed  in  Section  1.2,  L2  is  the  primary  frequency  band  of  interest. 
Therefore,  only  the  L2  band,  transmitting  the  new  civil  signal  (L2C),  the  legacy 
P(Y)-code  and  the  new  M-code,  is  considered  for  this  effort.  Furthermore,  given 
that  the  current  GPS  signal  power  levels  do  not  interfere  with  radar  systems  today, 
a  majority  of  the  analysis  is  focused  on  the  M-code. 

To  bound  the  simulation,  the  radar  system  is  modeled  as  a  single  pulse  radar 
and  the  target  is  assumed  stationary.  The  goal  of  the  simulation  is  to  determine  the 
effects  of  the  M-code  on  the  radar  system’s  ability  to  detect  a  target.  In  general, 
this  research  will  involve  a  worst-case  analysis  of  the  effect  of  GPS  M-code  on  an  Air 
Traffic  Control  (ATC)  radar. 

1.5  Thesis  Organization 

Chapter  2  presents  additional  information  on  the  GPS  M-code  signal  and  back¬ 
ground  information  on  radar  theory.  Chapter  3  presents  the  computer  simulation 
model,  simulation  results  and  analysis.  Finally,  Chapter  4  summarizes  the  results  of 
the  thesis,  discusses  limitations,  and  provides  recommendations  for  future  research. 
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II.  GPS  Signal  Structure  and  Radar  Theory 


2. 1  Overview 

This  chapter  presents  detailed  information  on  the  military  code  (M-code)  sig¬ 
nal  structure  and  a  discussion  on  radar  theory.  The  M-code  discussion  focuses  on 
signal  generation  as  compared  with  the  C/A  and  P-codes,  although,  they  will  not 
be  discussed  in  depth.  The  radar  theory  focuses  on  radar  pulse  characteristics,  the 
range  equation,  and  target  detection.  The  M-code  and  radar  theory  information 
presented  will  be  used  for  the  simulation  presented  in  Chapter  3.  For  a  good  review 
of  the  C/A  and  P-codes,  the  reader  is  directed  to  references  [7]  and  [10]. 

2.2  GPS  M-code 

Understanding  the  M-code  signal  structure  requires  some  insight  into  the  GPS 
upgrade  design  goals.  Some  of  these  goals  include  [4]: 

•  Better  jamming  resistance  than  the  current  P(Y)-code  signal  accomplished 
through  higher  transmit  power  without  causing  C/A-code  or  P(Y)-code  re¬ 
ceiver  interference. 

•  Compatibility  with  prevention  jamming  against  enemy  GPS  use. 

•  More  robust  signal  acquisition. 

•  Comparable,  if  not  better,  performance  to  the  current  P(Y)-code  signal. 

•  Coexistence  with  current  signals  on  LI  and  L2,  not  interfering  with  current  or 
future  military  user  equipment. 

•  Simple  and  low-risk  implementation  on  both  space  vehicles  and  in  future  user 
equipment  -  must  be  as  power  efficient  as  possible. 

To  achieve  these  goals,  the  M-code  was  designed  to  be  modulated  with  a  Bi¬ 
nary  Offset  Carrier  (BOC)  signal  having  a  subcarrier  frequency  of  10.23  MHz  and 
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spreading  code  rate  of  5.115  M  spreading  bits  per  second,  denoted  BOC(10. 23,5. 115) 
(abbreviated  as  BOC(10,5)).  BOC  modulation  details  are  sufficiently  covered  in  [5]. 

Mathematically  the  transmitted  M-code  signal  can  be  written  as  [8] 


SM(t)  =  y/2. P^dM(t)SW(t)PN5(t)cos(uL1,2t  +  6)  (2.1) 

where 


Pm 

=  Transmitted  M-code  power 

dM(t) 

=  M-code  data  stream  (25  or  100  bps) 

SW(t) 

=  10.23  MHz  Square  wave  carrier 

PN5(t) 

=  5.115  MHz  Pseudo  random  code 

UL1,2 

=  LI  or  L2  angular  frequency 

e 

=  Phase 

Figure  2.1  shows  an  example  of  SW(t),  PN5(t),  and  the  resultant  baseband 
M-code  sequence  SW(t)PN5(t). 


Time  (sec)  x1Q-s 

Figure  2.1  Example  of  the  M-code  PN  Sequence 


2-2 


The  main  result  of  the  BOC(10,5)  modulation  is  that  the  majority  of  the  power 
in  the  L2  frequency  band  is  concentrated  around  offset  frequencies.  These  frequencies 
are  located  10.23  MHz  above  and  below  the  carrier  frequency  as  shown  in  Figure  2.2, 
where  the  power  spectral  density  (PSD)  is  given  by  [14] 


PSDboc ~  fc 


(2.2) 


where  fs  =  10.23xl06  and  fc  =  5.115xl06. 


Figure  2.2  PSD  of  BOC(10,5)  M-Code  Modulation 


Figure  2.3  illustrates  a  comparison  of  the  M-code  spectrum  with  that  of  the 
C/A  and  P(Y)-code  signals,  where  all  signals  have  1  W  power.  From  Figure  2.3,  the 
M-code  modulation’s  spectrum  appears  to  be  distinctly  separated  from  that  of  the 
C/A  and  P(Y)-code  signals.  As  a  result,  the  M-code  may  be  received  at  high  power 
levels  without  degrading  the  performance  of  C/A  and  P(Y)-code  receivers — a  key 
design  goal  [6].  Table  2.1  shows  the  minimum  and  maximum  user- received  RF  signal 


2-3 


power  levels  for  the  M-code,  listed  by  satellite  production  version  [14].  As  stated 
previously,  the  effect  of  the  potential  increase  in  the  maximum  received  RF  power 
on  ATC  radar  systems’  detection  capability  is  the  primary  focus  of  this  thesis. 


Figure  2.3  PSD  Comparison  of  C/A,  P(Y),  and  M-code 


Table  2.1  Received  RF  M-code  Signal  Strength  Minima  and  Maxima 
(dBW)  [14] 


Production  Version 

Min 

Max 

Block  IIF 

-160 

-153 

Block  IIR-M 

-160 

-153 

Future  SVs 

-158 

-131 

2.3  Radar  Theory 

The  term  radar  is  an  acronym  for  RAdio  .Detection  And  Ranging  [17].  Radars 
operate  by  radiating  electromagnetic  energy  and  examining  the  reflected  energy  for 
detection  and  range  determination.  The  basic  radar  principle  is  illustrated  in  Fig¬ 
ure  2.4. 
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Figure  2.4  Illustration  of  the  Basic  Radar  Principle 


2.3.1  Radar  Pulse.  Most  objects  such  as  aircraft,  ships,  vehicles,  buildings, 
or  features  of  the  terrain  reflect  radio  waves,  much  as  they  do  light.  The  only 
difference  between  radio  waves  and  light  are  the  transmission  frequencies,  where  light 
is  much  higher.  The  reflected  energy  is  scattered  in  many  directions,  but  generally 
a  detectable  portion  of  it  is  scattered  back  to  the  original  transmission  point  [18]. 

In  most  radar  systems,  the  transmitter  and  receiver  share  a  common  antenna, 
as  shown  in  Figure  2.4.  In  this  configuration,  the  transmitter  may  interfere  with 
the  reception.  To  avoid  this  problem,  the  radio  waves  are  generally  transmitted  in 
pulses,  at  which  time  the  receiver  is  turned  off. 

The  following  terms  are  used  to  describe  the  pulse  transmission  as  illustrated 
in  Figure  2.5  [18]: 

•  Pulse  Duration  (r)  -  Length  of  the  transmitted  pulse 

•  Interpulse  Period  (T)  -  The  time  between  transmitted  pulses 

•  Pulse  Repetition  Frequency  (PRF)  -  Rate  at  which  pulses  are  transmitted 

(fr  =  y) 

During  each  interpulse  period,  T,  the  radar  radiates  energy  for  r  seconds  and 
“listens”  for  target  returns  for  the  rest  of  the  interpulse  period.  In  general,  the 
transmitted  signal  takes  the  form  of  a  pulsed  sine  wave  with  an  angular  frequency 
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Figure  2.5  Illustration  of  Pulse  Characteristics 


u0,  where  uia  is  the  carrier  frequency  or  operating  frequency.  For  long  range  military 
and  air  traffic  control  search  operations,  the  frequency  band  utilized  is  between  1.0- 
2.0  GHz,  classified  as  L-band. 

The  received  energy  returned  from  a  target  located  at  distance  R  will  have  a 
time  delay  given  by 

2  R 

A  t=  —  2.3 

c 

where  c  is  the  speed  of  light.  Assuming  the  reflected  energy  from  a  particular  pulse 
is  received  within  the  interpulse  period  T,  there  will  be  no  range  ambiguity,  and 
the  maximum  unambiguous  range  as  a  function  of  the  interpulse  period  or  pulse 
repetition  frequency,  fr,  is  given  by 


R 


max 


cT 

~2~ 


c 

2 Jr 


(2.4) 


As  an  example,  consider  the  case  shown  in  Figure  2.6.  Echo  1  represents  the  radar 
return  from  a  target  at  range  Ri  =  cAt/2  due  to  pulse  1.  Echo  2  could  be  interpreted 
as  the  return  from  the  same  target  dne  to  pulse  2,  or  it  may  be  the  return  from  a 
target  further  away  at  range  R2  dne  to  pulse  1.  The  fact  that  echo  2  could  be  from 
either  pulse  indicates  there  is  a  range  ambiguity  associated  with  echo  2  [12]. 
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Figure  2.6  Range  Ambiguity  Illustration 


2.3.2  Radar  Equation.  Following  the  development  of  [9],  the  radar  equation 
is  derived.  The  detection  range  of  a  radar  system  is  primarily  a  function  of  three 
parameters:  1)  transmitted  power,  2)  antenna  gain,  and  3)  receiver  sensitivity.  An 
increase  in  transmit  power  will  increase  the  radiated  energy,  which  will  result  in 
a  stronger  return  signal.  If  P  is  the  power  radiated  by  an  antenna  and  power  is 
radiated  uniformly  in  all  directions,  the  power  per  unit  area  or  power  density  at  a 
range  R  is  given  by 


Power  density 


P 

4nR2 


(2.5) 


where  AttR2  represents  the  area  of  a  sphere  with  radius  R. 


Most  radars  have  a  directive  antenna  that  radiates  power  in  a  specific  direction. 
This  directivity  is  accounted  for  by  incorporating  the  antenna  gain  into  Eq.  2.5.  The 
antenna  gain  at  a  particular  angle,  9,  is  defined  as  the  ratio  of  the  radiated  intensity 
at  6  to  the  radiation  intensity  of  a  uniformly  radiating  antenna.  The  angle  6  is 
illustrated  in  Figure  2.7  reproduced  from  [9].  At  6  —  0,  the  maximum  gain  of  an 
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antenna  is  related  to  its  physical  area  A  and  is  given  by 


_  4tt Apa 


(2.6) 


where  A  is  the  area,  pa  is  the  antenna  efficiency,  and  A  is  the  wavelength  of  the 
transmitted  wave.  Using  this  equation  for  antenna  gain,  the  power  density  equation 
given  in  Eq.  2.5  is  modified  to  provide  the  radiated  power  density  of  a  directive 
antenna  (Pda)  as  shown  by 


Pda  = 


PG 

AtcR2 


(2.7) 


BORESIGHT  AXIS 


Figure  2.7  Antenna  Gain  Characteristics 

When  the  radiated  energy  encounters  a  target,  a  portion  of  the  energy  is  re¬ 
flected  back  to  the  receiving  antenna.  The  amount  returned  is  dependent  on  the 
target  scattering  characteristics  or  the  Radar  Cross  Section  (RCS)  denoted  by  cr, 
which  has  units  of  area  (ft2,  m2).  The  RCS  takes  into  account  target  size,  orien¬ 
tation,  physical  shape,  and  material,  and  is  normally  measured  experimentally  [18]. 
Thus,  the  reflected  power  from  the  target  becomes 


Reflected  power  from  the  target 


PGa 
4t tR2 


(2.8) 
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Viewing  the  reflected  power  as  a  uniform  source  of  power  located  at  the  target, 


the  power  density  at  the  radar,  a  distance  R  away,  can  be  calculated  by 


Power  density  of  reflected  target  power  at  the  radar 


PGa 

(4vrP2)2 


(2.9) 


At  the  radar,  it  is  assumed  the  area  capturing  the  energy  is  equal  to  the  area 
of  the  receiving  antenna,  A.  Assuming  pa  =  1  for  math  simplicity,  the  power  of  the 
reflected  signal  at  the  radar  is  given  by 


Power  of  reflected  signal  at  the  radar 


PGaA 

(4vrP2)2 


(2.10) 


Solving  Eq.  2.6  for  A  and  substituting  in  Eq.  2.10  results  in  the  power  re¬ 
ceived  (Pr)  at  the  radar  as  shown  in  Eq.  2.11, 


PG2\2a 
(4vr  )3P4 


(2.11) 


The  target  return  power  is  only  half  of  the  story.  In  all  practical  applications, 
the  returned  signal  received  by  the  radar  is  corrupted  with  noise.  This  noise  can 
originate  within  the  receiver  itself  or  it  may  enter  the  receiver  through  the  receiving 
antenna.  The  M-code  signal  received  by  the  antenna  is  one  example  of  interference 
contributing  to  noise.  One  part  of  the  noise  generated  within  the  receiver  is  due  to 
the  thermal  agitation  of  electrons  in  the  receiver.  This  thermal  noise  has  the  value 
given  by 

Thermal  noise  =  kTBn  (2-12) 

where  k  =  1.38x10 ~2i joule/ degree  Kelvin  is  Boltzmann’s  constant,  T  is  the  noise 
temperature  in  degrees  Kelvin,  and  Bn  is  the  noise  bandwidth  in  Hz.  As  previously 
mentioned,  thermal  noise  is  only  one  part  of  the  receiver  noise.  To  account  for  all  of 
the  receiver  noise,  Eq.  2.12  is  multiplied  by  a  figure  of  merit  called  the  noise  figure, 
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Fn.  This  multiplication  results  in  the  total  receiver  noise  power 


N  =  FnkTBn 


(2.13) 


Using  the  noise  power  in  Eq.  2.13  and  denoting  the  target  power  received  by 
the  radar  in  Eq.  2.11  as  S,  the  signal-to-noise  ratio,  S/N,  forms  the  radar  equation 
given  as 


S  _  PG2X2a 
N  ~  (47r)3R4FnkTBn 


(2.14) 


where  S/N  is  generally  expressed  in  decibels  (dB). 

Since  the  target  return  signal  is  always  accompanied  by  noise,  the  receiver 
must  distinguish  the  signal  from  the  noise.  This  procedure  is  described  in  the  next 
section. 


2.3.3  Target  Detection.  The  following  discussion  on  target  detection  is 
adapted  from  [11].  All  of  the  material  presented  may  not  be  applicable  to  this 
research  effort;  however,  the  validity  of  this  analysis  will  be  addressed  at  the  end  of 
this  section. 

The  purpose  of  the  radar  receiver  is  to  process  the  signal  and  detect  the  ex¬ 
istence  of  a  target  in  the  presence  of  noise,  which  includes  the  interference  from 
the  M-code  signal  for  this  effort.  The  process  is  accomplished  by  the  detection  cir¬ 
cuit  within  a  receiver.  A  block  diagram  of  a  simple  detection  circuit  is  shown  in 
Figure  2.8  [11]. 


Bandpass 

Envelope 

©. 

Threshold 

Filter 

Detector 

Stage 

Figure  2.8  Simple  Detection  Circuit 


The  elements  that  make  up  the  circuit  include  a  bandpass  filter,  usually  at 
the  Intermediate  Frequency  (IF),  an  envelope  detector,  and  a  threshold  circuit.  The 
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threshold  circuit  compares  the  output  of  the  envelope  detector  with  a  predetermined 
threshold.  Whenever  the  envelope  exceeds  the  threshold,  a  target  is  assumed  to  be 
present.  This  concept  is  illustrated  in  Figure  2.9. 


Figure  2.9  Detection  Concept  [11] 

Figure  2.9  also  shows  one  noise  peak  that  would  be  mistaken  for  a  target  (the 
right-most  peak).  This  threshold  crossing  is  known  as  a  false  alarm.  In  addition, 
one  target  would  not  be  detected  because  it  is  below  the  threshold,  denoted  a  missed 
detection.  A  missed  detection  would  indicate  that  the  probability  of  detection  (Pd) 
is  less  than  one.  Lowering  the  threshold  would  capture  the  missed  target,  thereby 
increasing  the  detection  probability;  however,  it  may  also  increase  the  false  alarm 
probability  ( Pfa )•  One  way  to  increase  Pd  without  increasing  Pfa  would  be  to 
increase  the  signal  so  that  the  peak  would  be  higher  than  the  threshold,  i.e.,  increase 
the  signal-to-noise  ratio,  SNR. 

To  illustrate  the  discussion  of  Pd  and  Pfa-,  an  analysis  using  the  circuit  of 
Figure  2.8  will  be  described  for  the  detection  of  a  single-pulse.  Assuming  the  input 
consists  of  a  sinusoidal  signal  with  additive  white  Gaussian  noise  (AWGN),  two 
probability  density  functions  (PDF)  of  the  envelope  are  required:  a  noise  only  case 
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and  a  case  when  both  signal  and  noise  are  present.  Using  the  PDFs  and  selecting  a 
threshold  value  will  determine  the  PB  and  the  Pfa- 

Determination  of  these  two  joint  PDFs  requires  first  examining  the  noise. 
When  white  Gaussian  noise  is  passed  through  a  narrow  bandpass  filter,  the  out¬ 
put  noise  can  be  described  by 

n0(t )  =  X(t)cos{uct )  +  Y  (t)sin(uct)  (2-15) 

where  ujc  is  the  center  frequency  of  the  filter.  X(t)  and  Y(t)  are  two  independent 
random  variables  having  a  Gaussian  PDF  with  zero  mean  and  each  having  the  same 
variance  as  n0(t).  If  the  two-sided  noise  power  spectral  density  (PSD)  is  Na/ 2,  and 
the  bandpass  filter  has  a  rectangular  response  with  bandwidth  fB,  then  the  total 
average  noise  power,  Pn,  is  determined  by 

Pn  =  XHi)  =  THt )  =  rifc)  =  NJb  (2.16) 

where  the  bar  over  X,  Y,  and  n0  indicate  an  expected  value. 

The  signal,  which  is  a  sine  wave  of  duration  r,  and  frequency  cnc,  will  pass 
through  the  filter  virtually  unchanged  (assuming  fB  1/r).  At  the  bandpass 
filter  output  (point  A  in  Figure  2.8),  during  the  pulse  length  r,  the  signal  could  be 
described  by 

s0(t)  =  A  cos[u>ct  —  4>s\  =  a  cos(u!ct )  +  b  sin{uct )  (2-17) 

where 


A  =  (a2  +  b2)^2 


(2.18) 


and 


4>s  =  arctan— 
a 


(2.19) 
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The  combined  signal  plus  noise  at  the  output  of  the  filter  will  be  the  sum  of 


Eq.  2.15  and  Eq.  2.17 


e0(t) 

=  s0(t)+n0(t) 

=  [a  +  X  (t)]cos(uct]  +  [b  +  Y{t)]sin{ujct) 

=  r(t)cos[ojct  +  (pit)} 

(2.20) 

where 

r(f)  =  {[a  +  A'(t)]2  +  [6  +  V(f)]2}1/2 

(2.21) 

and 

b  +  Y(t) 

(pit)  —  arctan - — — 

w  a  +  X(t) 

(2.22) 

Referring  back  to  Figure  2.8,  a  linear  envelope  detector  will  yield,  at  its  output 

(point  B),  the  envelope 

r(t),  which  can  be  rewritten  as 

r(t)  =  [XKt)  +  U2(*)]1/2 

(2.23) 

where 

Xi  =  a  +  X(t),  Yi  =  b  +  Y(t) 

(2.24) 

Since  both  X  (t)  and  Y (t)  are  independent  zero-mean  Gaussian  random  vari¬ 
ables,  X\  and  Y\  are  Gaussian  random  variables  with  means  a  and  b  respectively. 
Both  are  also  independent  and  can  be  described  by  their  corresponding  PDFs  shown 
in  Eqs.  2.25  and  2.26. 


Pi(AG) 


^)V^exp 


-{Xi  -  a )2~ 
2(32 


(2.25) 


P2{Y1) 


^GeXP 


(Yi-b)2 1 
2  (32 


(2.26) 
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Both  X\  and  Y\  have  the  same  standard  deviation  (3,  which  is  equal  to  the  root 
mean  square  (RMS)  value  of  the  noise  n0(t)  shown  in  Eq.  2.27. 


P  =  K(t)]l/2  =  (nJb)1'2 


(2.27) 


Since  X\  and  Y\  are  independent,  their  two-dimensional  joint  PDF  is  equal  to 
the  product  of  Eqs.  2.25  and  2.26  and  is  given  by 


p(X1,Y1)  =  pi(Xi)p2(Yi) 


1  pxn  [ (X^ap+CYx-bp 
27T/32  exP  [  2/32 


(2.28) 


Now  a  transformation  to  two  other  variables  r  and  0  is  made  by  using  the 
relations 

X\  =  rcoscj)  Yl  =  rsincf)  (2.29) 


r  =  (X2  +  Y2)1^2,  0  =  arctan — 

X\ 


(2.30) 


The  joint  PDF  of  r  and  0  is  obtained  from  the  joint  PDF  of  X\  and  Y\  using 
the  transformation 

M\  _  (O  Ol) 


p(r  0)  =  17 

P[  |./(X1,Y'1)| 

where  J(X\,  Y\)  is  the  Jacobian  of  the  transformation 


dr 

dr 

Xi 

n 

dX1 

dY i 

r 

r 

d(j> 

dX, 

dcf> 

dY\ 

-Yl 

r2 

Xl 

r2 

T(V  V)  _  UI1  _  '  '  _ 

—  q,  ~  ,,  „  —  _ 

drp  dtp  —  \  1  X\  y 

dXi  dYi  r2  r2 

Substituting  Eqs.  2.28,  2.29,  and  2.32  into  Eq.  2.31  yields 


r  —  (r2  +  a2  +  b2  —  2ra  coscft  —  2rb  sirup) 

=  YpP  CXP  - 2/P - 


(2.31) 


(2.32) 


(2.33) 
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The  PDF  of  the  enveloper  can  be  obtained  by  integrating  p(r,  <fi)  over  all  the 


phases 


r»27r 


p(r )  =  /  p(r,  4>)  dcf) 
Jo 


(2.34) 


which  can  be  written  as 


r  —  r‘ 


p(r,  0)  =  exP 


+  A2) 


2ir(d2  ^  2p2 


/*27T 

rA 

/  exp 

Jo 

— cos(0  -  <ps) 

d<t>,  (2.35) 


where  A  and  <ps  are  the  amplitude  and  phase  of  the  signal  as  defined  in  Eq.  2.17. 
The  solution  of  the  integral 


1 

27T 


/*27 r 

rA 

/  exp 

Jo 

—  cos(0  -  0s) 

d0 


(2.36) 


is  known  as  the  modified  Bessel  function  of  order  zero,  Iq(x),  which  simplifies  the 


PDF  to 


T 

V (r)  =  jp  exp 


—  (r2  +  A2) 

.  W2  . 


(2.37) 


The  PDF  in  Eq.  2.37  is  called  Rician ,  named  after  S.O.  Rice.  In  developing 
Eq.  2.37  there  were  no  restrictions  on  the  signal  amplitude  A.  Therefore,  assuming 
A  =  0,  i.e. ,  no  signal,  will  yield  the  PDF  of  the  envelope  of  narrow-band  noise.  Given 
/o(0)  =  1  results  in  a  PDF  given  as 


Po(r)  =p(r)  U=o  =  ^ 


(2.38) 


which  is  the  Rayleigh  PDF.  Figure  2.10  shows  the  two  pdf  plots  given  by  Eqs.  2.37 
and  2.38.  The  signal  plus  noise  pdf  assumes  a  SNR  of  A2/{2j32)  =  8.  In  addition, 
a  threshold  level  Vr  is  shown  to  indicate  that  whenever  r  >  Vr ,  the  detection  of  a 
target  is  declared. 
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VT/P  r/p 


Figure  2.10  PDF  of  Detected  Noise  and  Signal  Plus  Noise  [11] 

To  determine  Pfa,  the  area  underneath  p0(r )  given  by  Eq.  2.38  needs  to  be 
evaluated  from  the  threshold  to  infinity  as  shown  in  Eq.  2.39. 

r  [  —  r  1  r_y2i 

Pfa  =  JvT  W  exp  M  Ar  =  exp  [W\  (239) 


Similarly,  to  determine  PD,  the  area  underneath  p(r)  given  by  Eq.  2.37  needs 
to  be  evaluated  from  the  threshold  to  infinity.  However,  the  integral  of  Eq.  2.37  does 
not  have  a  closed  form  solution  and  can  only  be  expressed  by  tabulated  functions. 
Yet,  if  the  SNR  is  high,  an  approximation  can  be  made.  Given 


n.rr,  S 

snr=n  =  W>>1 


(2.40) 


the  PDF  will  have  a  sharp  peak  near  r  =  A,  which  means  the  argument  of  the 
modified  Bessel  function,  rA,  can  be  replaced  with  A2.  Using  the  approximation 


I0(x )  «  (27ra:)  ^2exp(a:),  lCi 


(2.41) 
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yields,  for  large  SNR, 


p(r ) 


(2v r)1/2/? 


exp 


■0 r-Ay 

2(32 


(2.42) 


which  is  the  Gaussian  PDF.  The  integral  of  the  Gaussian  PDF  is  better  known, 
yielding 


Pd  =  P(Vt  <  r  <  oo) 


"1-er/(^-(SiVfl)1/2) 


(2.43) 


where  erf  is  the  error  function  defined  as 


erf(x )  =  /  e  t  dt 


7 T 


(2.44) 


While  the  target  detection  process  presented  in  this  chapter  utilizes  an  envelope 
detector,  the  simulation  that  will  be  presented  in  Chapter  3  employs  a  matched  filter 
detection  process.  The  information  on  the  envelope  detector  was  provided  to  show 
the  general  relationship  between  Pfa,  Pd ,  and  threshold  selection,  which  will  be 
applicable  in  the  matched  filter  implementation. 


2-4  Summary 

In  this  chapter,  the  mathematical  derivation  of  the  M-code  signal  was  pre¬ 
sented.  In  addition,  an  overview  of  radars  was  provided,  with  an  emphasis  on  de¬ 
tection,  particularly  how  the  probability  of  detection,  probability  of  false  alarm,  and 
SNR  interact.  This  working  knowledge  provides  a  conceptual  basis  for  the  simulation 
and  results  presented  in  Chapter  3. 
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III.  Simulation  Methodology,  Results,  and  Analysis 


3. 1  Overview 

This  chapter  presents  the  simulation  methodology  and  the  simulation  results 
of  M-code’s  effects  on  radar  detection.  A  simulation  description  is  presented  first, 
listing  the  assumptions  made  for  the  simulation.  Then  the  simulation  process  is 
described.  The  simulation  and  analysis  are  done  in  two  phases:  1)  the  baseline 
case  where  the  M-code  is  not  present  and  2)  the  interference  case  where  the  M-code 
(interferer)  is  present.  The  analysis  uses  detection  probability  (Pd)  curves  which 
show  the  detection  probability  versus  input  signal-to-noise  ratio. 

3.2  Simulation  Description 

As  stated  in  Section  1.2,  the  overall  objective  of  this  thesis  is  to  illustrate  the 
effects  the  M-code  signal  has  on  radar  systems.  In  particular,  this  simulation  focuses 
on  the  effects  the  M-code  signal  has  on  a  radar  system’s  detection  capability.  The 
simulation  is  implemented  in  two  phases.  The  first  phase  focuses  on  simulating  the 
detection  process  in  the  presence  of  thermal  noise  only,  without  the  presence  of  the 
M-code  signal.  The  second  phase  adds  the  M-code  signal  as  an  interference  source 
which  the  radar  system  must  overcome  in  its  detection  process. 

An  important  aspect  of  the  simulation  is  that  it  is  designed  to  look  at  a  worst 
case  radar  detection  scenario.  As  with  any  simulation,  assumptions  are  made  to  fa¬ 
cilitate  the  implementation.  The  assumptions  made  for  this  simulation  are  captured 
in  Figure  3.1  and  are  summarized  as  follows: 

1.  The  return  energy  is  from  a  stationary  airborne  point  target.  Doppler  and 

clutter  effects  are  not  considered. 
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Figure  3.1  Simulation  Scenario 


2.  The  simulation  is  done  for  a  single  pulse  return.  This  scenario  is  a  worst  case 
for  detection,  since  benefits  of  pulse  integration  are  not  included. 

3.  The  RF  GPS  M-code  signal  is  received  in-line  with  the  target  return  pulse. 
Physically,  this  implies  the  target  and  transmitting  GPS  satellite  are  in  direct 
line-of-sight  with  the  radar  (as  shown  in  Figure  3.1).  When  aligned  in  such 
a  manner,  the  received  M-code  signal  has  the  maximum  radar  antenna  gain 
applied  to  it  (so  it  is  worst-case). 

4.  The  radar  carrier  frequency  is  centered  at  one  of  the  main  lobes  of  the  M-code 
signal  (see  Figure  3.2).  This  location  simulates  a  maximum  received  M-code 
signal  power  condition  and  again  represents  the  worst-case  scenario. 

5.  Noise  is  modeled  as  additive  white  Gaussian  noise  (AWGN). 

6.  The  detection  threshold  value  is  set  using  only  AWGN  with  no  interference 
present. 

7.  The  M-code  PN  spreading  sequence  is  assumed  to  be  well  modeled  using  a  ran¬ 
dom  binary  sequence.  While  the  actual  M-code  is  classified,  this  is  a  reasonable 
assumption. 
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Figure  3.2  Radar  RF  Frequency  Chosen  to  Maximize  Received  M-code 
Signal  Power 

The  radar  system  model  for  this  simulation  is  based  on  the  Air  Route  Surveil¬ 
lance  Radar  Mode  4  (ARSR-4).  The  following  discussion  on  the  ARSR-4  system  is 
taken  from  Instruction  Book,  Field  Maintenance,  ARSR-4  System,  Type  FA- 10331 
Sections  1-10  [1]. 

The  ARSR-4  system  is  a  three-dimensional  long  range  radar  jointly  used  by  the 
Federal  Aviation  Administration  (FAA)  and  the  United  States  Air  Force  (USAF). 
The  system  provides  360-degree  azimuth  coverage  for  ranges  out  to  250  nautical 
miles,  at  a  height  of  up  to  100,000  feet,  and  for  elevation  angles  of  -7  to  +30  degrees. 

The  ARSR-4  radar  is  required  to  detect  very  small  targets  (a  =  0.1  m2  targets) 
at  ranges  up  to  92  nmi  and  a  =  2.2  m2  targets  at  ranges  up  to  200  nmi.  The  radar 
has  a  range  resolution  of  1/16  nmi  and  a  minimum  range  requirement  of  5  nmi.  To 
meet  these  requirements,  the  ARSR-4  radar  transmits  a  wide  pulse  (150  psecs)  made 
up  of  two  subpulses:  one  90  psecs  wide  and  the  other  60  psecs  wide.  The  two  pulses 
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are  transmitted  at  two  different  frequencies  separated  by  82.8572  MHz.  The  60  /zsec 
pulse,  transmitted  at  the  higher  frequency,  is  narrow  enough  to  meet  the  5  nmi  range 
requirement.  The  90  /zsec  pulse,  transmitted  at  the  lower  frequency,  provides  the 
additional  energy  needed  to  meet  the  larger  range  detection  requirements.  To  meet 
the  range  resolution,  the  return  echo  is  compressed  to  0.772  yusecs  which  implies 
that  each  transmit  subpulse  is  digitally  phase-coded  to  create  Non-linear  Frequency 
Modulation  (NLFM)  with  1.3  MHz  deviation.  The  radar  transmits  the  90  /jl sec 
pulse  first,  followed  by  the  60  p.sec  pulse.  The  operating  parameters  for  the  ARSR-4 
system  are  summarized  in  Table  3.1. 


Table  3.1  Operating  Parameters  for  the  ARSR-4  System 


Parameter 

Value 

Peak  Transmit  Power 
Average  Power 

Waveform  Duty  Cycle 
Antenna  Transmit  Gain 
Antenna  Receive  Gain 
Frequency 

PRF 

Scan  Time 

Azimuth  Beamwidth 

63.765  kW 

2.55  kW 

4.32  % 

37.7  dB 

40.91  dB,  max 

1215-1400  MHz,  cliplex  operation 
288  Hz 

12  secs 

1.4  degrees 

As  previously  stated,  in  this  simulation  the  ARSR-4  system  is  simplified  to 
be  a  single,  uncompressed-pulse  radar  system,  which  is  a  worst-case  approach.  Use 
of  the  NLFM  pulse  compression  should  improve  detection  over  the  results  shown  in 
this  research. 

3.3  Simulation  Process 

Figure  3.3  illustrates  the  radar  receiver  model  used  for  the  simulations.  The 
received  signal  (either  signal  plus  noise  or  signal  plus  noise  plus  interference)  is 
initially  filtered  by  an  RF  Band  Pass  Filter  (BPFrf).  The  filtered  signal  then  goes 
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through  a  two-step  down-conversion  process  to  take  the  signal  from  an  RF-level  down 
to  an  IF-level  suitable  for  detection  processing. 


V 


Figure  3.3  Receiver  Model  (Shaded  Area  Represents  Portion  Simulated 
in  this  Research) 

Given  the  assumption  the  radar  frequency  is  centered  at  a  main  lobe  of  the 
M-code  signal,  it  is  assumed  the  effects  of  the  down-conversion  process  equally  affect 
the  returned  radar  pulse  and  the  portion  of  the  M-code  signal  received  by  the  radar 
receiver.  Thus,  front-end  filtering  effects  are  neglected  and  the  resulting  SNR  after 
the  second  local  oscillator,  LO2,  is  the  variable  simulation  parameter.  (If  the  radar 
frequency  is  not  at  a  main  lobe  of  the  M-code  signal,  the  M-code  signal  would  have 
less  impact  than  in  the  worst-case  scenario  that  was  simulated). 

The  remaining  discussion  on  the  simulation  process  focuses  on  the  shaded  re¬ 
gion  of  Figure  3.3.  The  simulation  is  broken  up  into  two  cases-the  baseline  case 
(no  M-code  signal  present),  and  the  M-code  case  (with  M-code  present).  Figure  3.4 
presents  a  simulation  flow  diagram  which  describes  the  individual  steps  of  the  sim¬ 
ulation.  The  shaded  box  represents  the  baseline  case,  while  the  rest  of  the  diagram 
illustrates  the  steps  that  are  followed  for  the  interference  case. 

3.3.1  Baseline  Case.  The  block  diagram  for  the  baseline  case  of  the  sim¬ 
ulation  is  shown  in  Figure  3.5.  The  BPFip  inputs  are  the  return  pulse,  s(t),  and 
the  noise  term,  n(t).  The  return  pulse  is  generated  as  a  pulsed  sinusoid  at  the  IF 
frequency  with  duration  r,  while  the  noise  is  zero-mean,  additive  white  Gaussian 
noise  (AWGN). 
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Figure  3.4  Simulation  Flow  Diagram  to  Determine  the  Effects  of 
M-code  on  Radar  Detection.  Shaded  Region  Represents  the 
Baseline  Case. 


n(t) 

Figure  3.5  Baseline  Simulation  Block  Diagram 
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The  signal-to-noise  ratio  going  into  BPFif,  SNRiF,  is  determined  by  taking  the 
ratio  of  signal  power,  Ps,  and  noise  power,  Pn  calculated  by 

1  jv  „ 

Ps  =  E  s2(tj)  (3.1) 

j=l 

1  N  „ 

Pn  =  W7  S  n 2(tj)  (3.2) 

IV  7=1 

given  a  normalized  lfl  load.  N  is  the  number  of  samples  in  either  the  signal  or  the 
noise  vectors.  Since  the  simulation  uses  sampled  data,  the  Ps  and  Pn  calculated  by 
Eqs.  3.1  and  3.2  represent  average  signal  and  noise  power.  In  addition,  in  determining 
the  noise  power,  a  noise  bandwidth  given  by  1/At  is  implied,  where  At  is  the  time 
difference  between  two  adjacent  samples  of  the  signal  or  noise  vector.  For  this 
simulation,  At  is  1.25  x  10-7  secs,  which  gives  an  implied  noise  bandwidth  of  8  MHz. 

The  BPFip  bandwidth  for  this  simulation  was  set  to  1/r  =  1/90  /xsecs  or 
approximately  11  kHZ,  which  is  the  3.0  dB  (or  half-power)  bandwidth.  This  value 
is  typical  for  radar  receivers  of  this  type  [9].  The  filter  selected  for  this  simulation 
is  an  8^-order  Chebyshev  bandpass  filter  with  0.1  dB  bandpass  ripple.  The  filter 
responses  of  the  Chebyshev  filter  and  an  ideal  filter  are  shown  in  Figure  3.6. 

The  BPFip  filtering  operation  yields  both  filtered  signal  and  noise  components 
that  are  input  into  the  matched  filter.  The  matched  filter  is  the  heart  of  the  detection 
process.  The  matched  filter  input  SNR,  SNR,D,  is  determined  in  the  same  manner 
used  to  determine  SNRIF,  except  the  samples  are  now  taken  after  the  BPFif. 

The  matched  filter  is  used  to  “match”  the  incoming  signal  with  a  known  refer¬ 
ence  pulse.  For  this  simulation,  the  original  signal,  s(t),  is  modeled  as  the  reference 
pulse.  The  matched  filter  can  be  implemented  either  as  a  correlation  operation  or 
a  convolution  operation.  The  latter  is  chosen  for  this  effort.  Convolution  initially 
requires  reversing  either  the  incoming  signal  or  the  reference  signal  in  time.  The 
process  involves  convolving  the  incoming  (filtered)  signal  with  the  reference  signal 
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Figure  3.6  Filter  Responses  of  an  8^-Order  Chebyshev  Filter  and  an 
Ideal  Filter 

by  sliding  the  reference  signal  in  time  until  the  correlation  is  maximized.  Since  it 
is  assumed  the  time  delay  of  the  return  pulse  is  known  and  the  return  signal  and 
reference  signal  sequences  are  the  same  length,  the  correlation  is  maximized  when 
the  two  sequences  completely  overlap  at  time  t  —  r,  where  r  is  the  pulse  width. 
Thus,  the  maximum  correlation  will  occur  at  the  Nth  sample  out  of  the  matched 
filter,  since  the  signal  and  noise  are  both  N  samples  long. 

Since  the  Nth  sample  is  the  only  value  of  concern  (because  the  exact  value  of  the 
return  signal  delay  is  assumed  to  be  known),  the  time  consuming  convolution  process 
was  simplified  for  the  simulation  by  an  equivalent  one-step  process.  The  signal  and 
noise  were  multiplied  sample  by  sample  and  individual  sample  products  were  then 
added.  The  results  of  this  multiplication  and  summation  operation  provides  the 
same  correlation  value  as  the  full  convolution  method  described  above. 

The  matched  filter  output  is  then  compared  to  a  threshold  value  to  determine 
detection.  Any  value  above  the  threshold  indicates  target  detection,  while  a  value 
below  indicates  no  target  present.  As  presented  in  Chapter  2,  the  probability  of 


3-8 


detection  (Pd)  and  the  desired  probability  of  false  alarm  rate  (Pfa)  depend  on  the 
threshold  value.  The  desired  threshold  value  is  determined  by  fixing  a  value  for 
Pfa  and  passing  the  noise  samples  through  the  matched  filter.  To  be  statistically 
perfect,  an  infinite  number  of  noise  realizations  (Nneai)  are  required,  which  is  not 
feasible.  Therefore,  the  Pfa  is  used  to  determine  the  minimum  number  of  simulated 
realizations  needed  to  be  statistically  accurate,  following  the  rule-of-thumb 

NReal  —  75 —  (3.3) 

-UFA 

The  Nth  sample  for  each  noise  realization  passed  through  the  matched  filter 
is  stored  as  a  test  statistic  (Z),  which  is  used  in  establishing  a  detection  threshold. 
Once  all  the  noise  test  statistics  are  collected,  they  are  sorted  and  the  threshold 
value  is  assigned  to  the  Njteai  —  (Pfa  x  Nneai )  largest  Z  value.  For  example,  for  a 
Pfa  —  0.01,  1000  noise  realizations  are  required  so  that  the  990t/l  largest  Z  value  is 
assigned  as  the  threshold.  This  example  is  illustrated  in  Figure  3.7. 

Threshold  Level  Determination 

60  | - ( - j - 1 - 1 - 1 - | - 1 - 1 - 1 - 


-40  -  .  ; 

-60 - 1 - 1 - 1 - l! - 1 - 1 - 1 - 1 - 1 - 

0  100  200  300  400  500  600  700  800  900  1000 

Realization 

Figure  3.7  Threshold  Determination  Example 
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The  matched  filter  processing  is  then  repeated  with  a  signal  added  to  each  noise 
realization  to  produce  a  new  test  statistic,  Z' .  The  number  of  Z'  values  exceeding 
the  pre-determincd  threshold  value  are  counted  and  divided  by  the  total  number  of 
Z  values  to  estimate  Pd-  The  estimated  Pd  is  valid  for  a  given  Pda  and  SNRd. 
Therefore,  to  generate  a  complete  detection  picture  (Pd  versus  SNR  plot)  requires 
the  process  be  repeated  over  a  range  of  desired  SNRD  values,  where  each  SNRD 
produces  a  Pd  for  a  constant  Pfa- 

Typically,  radar  applications  require  a  Pfa  on  the  order  of  1CT6.  However,  due 
to  computational  resource  limitations  and  time  constraints,  all  subsequent  research 
results  (including  the  interference  case)  are  based  on  PFa  =  10~2. 

Figure  3.8  shows  a  detection  curve  for  the  baseline  simulation  (no  M-code 
present).  The  signal  used  to  generate  the  curve  is  a  pulsed  sine  wave  of  duration 
r  =  90  /r secs  at  an  IF  frequency  of  400  kHz,  consistent  with  ARSR-4  specifica¬ 
tions  [1].  The  noise  is  generated  as  zero-mean,  AWGN  with  an  unit  magnitude 
average  power.  The  value  of  SNRD  was  changed  by  varying  the  pulse  amplitude 
while  keeping  the  noise  power  constant.  The  curve  was  generated  by  taking  the 
noise  samples  and  passing  them  through  the  8f,l-order  Chebyshev  BPFif  and  the 
matched  filter  detection  process  to  establish  a  detection  threshold  for  a  Pfa  of  0.01. 
The  signal  amplitudes  were  then  adjusted  for  varying  SNRD  values,  and  the  sig¬ 
nal  plus  noise  was  then  passed  through  the  BPFif  and  matched  filter  to  determine 
detection. 

As  Figure  3.4  indicates,  operating  points  are  selected  from  the  baseline  simu¬ 
lation  to  serve  as  the  no  interference  reference  points.  From  Figure  3.8,  selecting  a 
value  for  Pd  fixes  a  value  for  SNRd-  Table  3.2  lists  five  operating  points  that  will 
be  used  as  part  of  the  analysis  in  the  next  section. 

The  generation  of  Figure  3.8  and  the  selection  of  operating  point  parameters 
completes  the  baseline  (i.e. ,  no  M-code  interference)  case.  The  next  section  presents 
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Figure  3.8  Probability  of  Detection  Curve  for  a  Pfa  =  0.01 


Table  3.2 


Baseline  Simulation  Operating  Points  Used  for  Analysis 


Pd 

SNRd  (dB) 

0.9 

6.06 

0.8 

4.83 

0.7 

3.81 

0.6 

2.84 

0.5 

1.86 
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the  second  phase  of  the  simulation,  which  includes  the  modeling  of  a  return  pulse 
with  the  addition  of  M-code  interference. 


3.3.2  Signal  Plus  Noise  Plus  M-code  Case.  The  block  diagram  for  the 
second  phase  of  the  simulation  is  shown  in  Figure  3.9.  The  only  difference  between 
Figure  3.9  and  Figure  3.5  is  the  addition  of  the  M-code  signal  term,  M(t). 


Figure  3.9  M-code  Interference  Simulation 


In  the  baseline  case,  the  noise  power  was  kept  constant  while  the  return  pulse 
amplitude  was  varied  to  obtain  SNRn.  For  the  interference  case,  the  return  pulse 
power  is  kept  constant  and  the  noise  samples  are  varied  in  magnitude  to  obtain  the 
desired  SNRd  value.  The  return  pulse  power  is  held  constant  because  it  is  important 
to  model  the  relative  power  levels  between  the  radar  and  M-code  signals  correctly. 
This  power  relationship  is  maintained  by  calculating  the  radar  return  signal  power 
and  the  M-code  power  at  the  receive  antenna  output.  (All  other  operations  occurring 
after  that  point  are  the  same  for  both  signals  and  will  not  significantly  affect  the 
relative  power  levels  between  the  radar  return  and  M-code  signal). 

The  radar  return  power  after  the  receive  antenna  is  calculated  by  Eq.  2.11. 
For  a  target  with  a  a  =  2.2  m2  cross-section  at  a  detection  range  of  200  nmi,  the 
received  power  was  calculated  with  the  parameters  listed  from  Table  3.1  (using  a 
carrier  frequency  of  1217.37  MHz)  (see  Figure  3.2).  The  target  cross  section  value 
and  the  detection  range  are  consistent  with  the  ARSR-4  system  operation  [1],  Since 
200  nmi  is  the  maximum  detection  range,  the  return  signal  power  is  a  minimum  for 
a  target  located  at  this  distance.  It  is  assumed  that  the  M-code  signal  will  have  the 
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most  interfering  affect  when  the  return  signal  power  is  at  a  minimum.  Thus,  the 
received  signal  power  after  the  antenna  is  calculated  as  follows 


PG2X2a 
(4? r)3R4 

(63,765  x  103W)(10^)(10^)(12gS4z)2(2.2m2) 

(4tt)3  (200nmi  •  ^g^f)4 
=  16.56  x  10~15W  («  -138  dBW) 


for  the  target  of  interest.  (Note  that  the  G 2  term  is  the  product  of  two  gains,  because 
there  are  different  gains  for  the  transmit  and  receive  antennas.) 

Since  the  return  power  is  fixed,  the  noise  sample  amplitudes  (power)  are  ad¬ 
justed  to  provide  the  variable  parameter  SNRIF.  For  every  value  of  SNRIF,  the 
original  unit-power  noise  samples  are  scaled  by  a  gain  factor  given  by 


Gain  Factor 


I  Pr 

1  /  SNRjp 

V  lO^cr 


(3.4) 


The  gain  factor  was  derived  by  solving  the  basic  SNRiF  equation  given  by 
Eq.  3.5  for  Pn  and  taking  the  square  root. 

SNRif  =  101oglo  fp-)  (3.5) 

where  Ps  and  Pr  are  equivalent  terms. 

By  properly  adjusting  the  noise  samples  for  the  given  return  pulse  to  the  correct 
SNRif,  the  SNRd  value,  after  the  signal  plus  noise  is  passed  through  the  BPFif,  will 
be  correct  as  well.  This  scenario  is  illustrated  in  Figure  3.10  which  shows  the  required 
SNRif  needed  to  get  a  corresponding  SNRD  for  a  given  PD. 

Since  the  noise  power  has  changed,  a  new  detection  threshold  value  is  deter¬ 
mined.  The  process  to  establish  a  new  detection  threshold  is  the  same  presented 
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Figure  3.10  Pd  versus  SNRd/SNRif  for  Pfa  —  0.01 


for  the  first  phase,  except  the  noise  in  this  case  is  the  scaled,  filtered  noise  signal. 
Once  the  return  pulse  signal  and  noise  are  correctly  set,  the  M-code  received  power 
is  varied  to  determine  at  what  received  power  level  the  M-code  begins  to  effect  the 
established  operating  point  Pd- 

The  M-code  signal  is  generated  using  Eq.  2.1.  As  previously  stated,  the  M-code 
PN  sequence  is  assumed  to  be  random  and  is  generated  as  such.  In  addition,  the 
M-code  data  stream,  dM,  is  assumed  to  be  constant  for  the  given  signal  duration. 
(This  assumption  is  reasonable  because  the  data  bit  interval  is  much  longer  than  the 
return  pulse  duration  r.) 

The  generated  baseband  M-code  sequence  is  used  as  the  received  M-code  signal. 
The  goal  is  to  simulate  a  filtered  1/r  bandwidth  portion  of  the  GPS  M-code,  where 
both  the  peak  of  the  M-code  signal  and  the  filter  are  centered  at  IF  (400  kHz). 
However,  since  the  bandwidth  of  the  M-code  signal  (24  MHz)  is  much  wider  than 
the  center  of  the  BPFjf  filter  (400  kHz),  a  pre-filtering  step  is  done  to  approximate 
the  filter  effects  that  would  normally  be  in  the  front  end  of  the  radar  receiver.  This 
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operation  eliminates  frequency  aliasing  of  the  M-code  signal  about  zero  frequency 
that  would  occur  if  the  signal  was  simply  formed  with  the  main  lobe  centered  at 
400  kHz.  First,  the  baseband  M-code  sequence  is  generated  at  an  IF  frequency  of 
20.23MHz.  The  resulting  IF  signal  is  passed  through  a  bandpass  filter  centered  at 
10MHz  with  a  bandwidth  of  500KHz  to  capture  the  lower  main  lobe  of  the  M-code 
signal.  Finally,  the  filtered  M-code  signal  is  down-converted  to  the  center  frequency 
of  the  IF  filter  located  at  400KHz.  It  is  assumed  that  this  procedure  is  equivalent  to 
most  downconversion  and  filtering  processes  that  are  found  in  radar  receiver  front 
ends  and  will  result  in  capturing  the  maximum  M-code  signal  power  as  previously 
shown  in  Figure  3.2.  Once  the  M-code  signal  is  passed  through  the  BPFif,  it  has 
the  desired  1/r  bandwidth. 

Initially,  the  received  M-code  power  is  scaled  by  the  maximum  receive  antenna 
gain  given  in  Table  3.1,  since  the  return  pulse  power  was  calculated  using  this  gain. 
The  scaled  M-code  signal  (interferer)  is  then  added  to  the  return  pulse  signal  and 
noise.  All  three  are  then  passed  through  the  BPFif  and  the  matched  filter  detection 
process  which  produces  the  Pd  of  interest. 

In  accordance  with  proposed  M-code  power  levels,  the  M-code  signal  was  added 
at  received  power  levels  ranging  from  -160  dBW  to  -130  dBW,  where  the  received 
power  is  considered  to  be  the  power  incident  on  the  antenna.  Applying  the  maximum 
receive  antenna  gain  of  40.91  dB,  the  M-code  signal  prior  to  the  BPFif  was  scaled 
from  -160  dBW  and  -130  dBW  up  to  -119.09  dBW  and  -90.09  dBW,  respectively. 

Figure  3.11  shows  several  Pd  curves  at  different  M-code  interference  values  for 
a  single  M-code  signal.  (The  effects  of  varying  the  phase  or  PN  code  of  the  M-code 
signal  will  be  described  later).  Note  that  the  x-axis  shows  the  signal-to-noise  ratio 
(SNRd)  not  the  signal-to-noise-plus-interference  ratio.  The  results  indicate  that 
there  is  a  minimal  effect  on  the  detection  probability  for  the  minimum  specified 
received  M-code  power  of  -160  dBW.  However,  as  this  particular  M-code  signal  is 
increased  to  its  maximum  value  of  -131  dBW,  a  larger  SNRd  is  required  to  obtain  the 
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same  Pd-  For  a  Pd  =  0.9,  the  increase  in  required  SNRd  going  from  an  interfering 
signal  power  of  -160  dBW  to  -130  dBW  is  about  7.44  dB.  An  alternative  way  to 
view  the  results  is  presented  in  Figure  3.12,  which  shows  Pd  as  a  function  of  received 
M-code  power  for  the  given  operating  points.  Again,  at  a  received  M-code  power  of 
-160  dBW,  there  are  minimal  effects,  while  at  the  maximum  specified  M-code  signal 
power,  the  PD  for  all  operating  points  dropped  considerably. 


Figure  3.11  Detection  Curves  for  Varying  Levels  of  Received  M-code 
Power  for  Pfa  =  0.01 


Since  the  radar  was  modeled  as  a  single  pulse  system,  the  increase  in  required 
SNRd  for  an  increase  in  received  M-code  signal  power  seems  appropriate.  Most 
radar  systems,  including  the  ARSR-4,  integrate  over  multiple  pulses,  which  should 
(in  theory)  decrease  SNRd  to  produce  the  same  Pd-  Additional  improvements  can 
be  obtained  by  a  radar  waveform  providing  orthogonality  to  the  BPSK  signal.  Im¬ 
plementing  NLFM  will  provide  orthogonality  and  is  a  topic  that  will  be  revisited  as 
a  recommendation  for  further  study. 
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Figure  3.12  Pd  vs  Received  M-code  Power  for  Pfa  =  0.01 

Thus  far,  the  results  assumed  the  M-code  was  received  perfectly  in-phase  with 
the  radar  return.  Figures  3.13  -  3.16  illustrates  the  effects  of  phase  offsets  ranging 
from  — 7T  to  7T  in  increments  of  22.5°  for  received  M-code  powers  of  -160,  -150,  - 
140,  and  -130  dBW.  Figure  3.17  shows  the  case  for  the  specified  -131  dBW  power 
level.  In  all  the  figures,  the  dashed  line  is  included  as  reference  and  represents  the  no 
interference  detection  curve.  For  received  power  of -160  dBW  (Figure  3.13),  changing 
phase  offset  values  had  little  effect  as  all  the  curves  are  tightly  bunched  together. 
However,  as  the  received  power  increases  up  to  -130  dBW  (Figure  3.16),  the  effects 
are  dramatic.  Depending  on  the  phase  offset,  the  curve  shifts  substantially,  so  much 
that  at  certain  offsets,  the  Pd  never  increases  above  a  nominal  value. 

The  curves  in  Figures  3.11  and  3.12  were  generated  using  the  same  randomly 
generated  M-code  sequence.  When  different  sets  of  random  sequences  were  applied, 
the  curves  shifted.  Examples  of  the  shift  are  illustrated  in  Figures  3.18  and  3.19 
which  show  plots  of  10  different  M-code  samples  for  received  M-code  power  levels  of 
-160  dBW  and  -131  dBW,  respectively.  Again,  the  dashed  line  is  included  as  reference 
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Figure  3.13  Effects  of  Phase  Offsets  on  PD  for  Received  M-code  Power 
of  -160  dBW  and  Pfa  —  0  01 


Figure  3.14  Effects  of  Phase  Offsets  on  Pd  for  Received  M-code  Power 
of  -150  dBW  and  Pfa  =  0.01 
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Figure 


Figure 


.15  Effects  of  Phase  Offsets  on  PD  for  Received  M-code  Power 
of  -140  dBW  and  PFa  =  0.01 


.16  Effects  of  Phase  Offsets  on  PF  for  Received  M-code  Power 
of  -130  dBW  and  PFa  =  0.01 
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Figure  3.17  Effects  of  Phase  Offsets  on  Pd  for  Maximum  Specified  Re¬ 
ceived  M-code  Power  of  -131  dBW  and  PFa  =  0.01 

and  represents  the  no  interference  detection  curve.  As  observed  with  changing  the 
phase  offset,  at  the  minimum  received  power  level  (Figure  3.18),  the  effects  of  using 
different  M-code  PN  sequences  were  minimal.  At  the  maximum  received  power  level 
(Figure  3.19),  the  effects  were  severe.  The  difference  in  required  SNRd  at  Pd  =  0.8 
spanned  25  dB.  Note  that  there  are  a  total  of  460  PN  values  over  the  90  /rsec  pulse 
duration.  Analyzing  the  sequences,  4  of  the  10  had  more  +1  values  than  -1  values, 
but  there  was  no  correlation  between  this  fact  and  the  generated  curves  1 .  Therefore, 
it  seems  that  in  addition  to  the  phase,  the  M-code  PN  sequence  over  the  90  /x sec  pulse 
duration  appears  to  be  an  important  factor  in  determining  the  effects  on  detection 
at  the  maximum  specified  received  M-code  power. 

For  both  the  phase  offset  and  random  sequence  results,  the  fact  that  the  curves 
can  shift  either  to  the  left  or  to  the  right  would  indicate  the  M-code  signal  can  either 
add  or  subtract  to  the  return  signal,  which  is  held  constant.  A  curve  shift  to  the 

1This  is  probably  due  to  the  fact  that,  because  the  signal  is  sampled,  there  is  no  guarantee  that 
each  +1  (or-1)  M-code  chip  will  increase  or  decrease  the  matched  filter  output  by  the  same  amount. 
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Figure  3.18  Comparison  of  PD  Curves  Using  Different  Random  M-code 
Sequences  for  Received  M-code  Power  of  -160  dBW 


Figure  3.19  Comparison  of  Pd  Curves  Using  Different  Random  M-code 
Sequences  for  Received  M-code  Power  of  -131  dBW 
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left  indicates  an  increase  in  false  alarms  since  the  M-code  signal  is  forcing  more 
threshold  crossings  by  effectively  increasing  the  return  signal  samples  (add  effect). 
A  curve  shift  to  right,  to  the  point  where  no  detection  is  ever  achieved,  indicates 
the  M-code  signal  is  reducing  the  return  signal  samples  so  that  the  Pd  is  reduced 
(subtract  effect).  By  changing  the  phase  or  using  a  different  random  sequence,  the 
M-code  signal  had  different  effects  (add  or  subtract)  to  the  return  signal  on  a  sample 
by  sample  basis,  which  led  to  the  results  presented  above. 

3-4  Summary 

In  this  chapter,  the  simulation  methodology,  simulation  results  and  analysis 
to  show  the  effects  of  the  M-code  on  radar  detection  capability  was  presented.  Two 
phases  were  discussed  1)  a  baseline  case  where  only  signal  and  noise  were  considered 
and  2)  an  interference  case  where  the  M-code  signal  was  added.  The  results  indi¬ 
cate  for  the  minimum  specified  M-code  power,  the  effects  on  detection  arc  minimal. 
However,  at  the  maximum  specified  M-code  power,  detection  is  greatly  affected. 
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IV.  Conclusion  and  Recommendations 


4.1  Conclusion 

This  thesis  presented  the  theory,  simulation  methodology,  and  results  of  simu¬ 
lating  the  detection  process  of  a  radar  system  and  the  effects  the  M-code  signal  has 
on  that  process.  The  simulation  looked  at  a  worst  case  scenario  with  the  following 
key  assumptions: 

•  The  ARSR-4  was  modeled  as  a  single,  uncompressed  pulse  radar  system  with 
no  integration  benefits. 

•  The  return  signal  is  from  a  stationary  target  located  at  the  maximum  target 
detection  range.  This  scenario  results  in  the  lowest  (worst-case)  return  power. 

•  The  target  and  the  GPS  satellite  broadcasting  the  M-code  signal  are  in  direct 
line-of-sight  with  the  radar.  The  maximum  radar  antenna  gain  is  applied  to 
the  M-code  in  this  configuration. 

•  The  detection  threshold  value  is  set  using  only  AWGN  with  no  interference 
present. 

Additional  assumptions  were  made  and  are  listed  in  Section  3.2. 

For  a  single  pulse  radar  system,  the  results  showed  that  M-code  effects  on 
the  radar’s  detection  capability  were  minimal  when  the  received  signal  power  was 
at  the  minimum  specified  value  of  —160  dBW.  For  the  maximum  specified  M-code 
power  of  —131  dBW,  the  results  indicate  that  the  M-code  signal  will  effect  this  radar 
system.  The  severity  of  the  effects  were  shown  to  be  dependent  on  the  selection  of 
the  random  M-code  sequence  as  well  as  phase  offsets.  In  the  simulation,  10  different 
random  sequences  were  applied  and  the  results  showed  a  difference  in  required  SNRd 
of  approximately  25  dB  for  cases  when  a  detection  curve  was  generated.  In  three 
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cases,  the  Pd  never  increased  above  a  nominal  value.  The  effects  of  phase  offsets  had 
similar  results.  Note  that  these  results  only  apply  to  a  worst-case  scenario  for  a  single 
pulse  as  simulated  for  this  research.  Further  modeling  of  the  ARSR-4  radar  system, 
to  include  use  of  NLFM  pulse  compression  and  pulse  integration,  would  be  necessary 
before  definitively  concluding  that  the  M-code  does  or  does  not  significantly  affect 
radar  system  performance. 

4-2  Recommendations 

As  stated,  the  worst  case  scenario  assumed  several  things  about  the  radar  sys¬ 
tem  and  the  M-code  signal.  Therefore,  to  augment  this  research  effort,  the  following 
list  of  recommendations  should  be  explored  for  a  follow-on  effort: 

1.  The  simulation  should  be  executed  using  a  Pfa  =  10~6.  As  discussed,  most 
systems,  including  the  ARSR-4,  require  this  Pfa-  Due  to  time  constraints,  this 
research  was  performed  using  a  Pfa  —  ICR2. 

2.  Model  the  radar  system  as  a  NLFM  pulse  integrating  radar  detecting  a  moving 
target  within  a  cluttered  environment.  As  briefly  mentioned  in  Chapter  3, 
providing  orthogonality  between  the  radar  pulse  and  the  M-code  in  the  receiver, 
i.e.  matched  filter,  will  have  a  significant  impact  on  the  M-code  interference 
rejection  characteristics  of  the  radar. 

3.  Obtain  samples  of  the  real  M-code  PN  sequence,  since  the  M-code  PN  sequence 
was  generated  as  a  random,  binary  sequence  for  the  simulation.  As  the  results 
indicated,  the  selection  of  the  random  sequence  influenced  how  much  effect  the 
M-code  had  on  the  detection  capability. 

4.  The  assumptions  made  for  the  radar  receiver  front  end  needs  to  be  further 
developed.  The  effects  of  the  two-step  downconversion  were  neglected  and 
should  be  included. 
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5.  Evaluate  other  radar  systems  operating  in  the  L2  frequency  band.  Other  sys¬ 
tems  will  have  operating  parameters  different  than  the  ARSR-4  system. 

6.  Figure  3.7  showed  an  example  of  determining  a  threshold  value  where  the  test 
statistics  were  both  positive  and  negative.  Typically,  test  statistics  should  not 
be  allowed  to  go  negative.  Therefore,  the  magnitude  or  squared  magnitude  of 
the  test  statistics  should  be  employed. 
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